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•  Given  reel  numbers  a2«...ta^^ w*-«M^:4tttere»ted^4if  the 


classic  problem  of  the  error  In  computing  S 


•R) 


when  the  sum 


I  ^  p  *\  /'*  \  '  "  _ 

Is  computed  by  ■  la.  where  (a.  I  Is  the  nearest  Integer  to 

_ X  y  X.  (.  i)'j( 

fa..'  We  shall  first  study  this  error  as  a  function  of  a  A  shift, 

l.e.,  when  all  numbers  a^  are  each  shifted  A  and  then  rounded; 


I  (•i-A) 
1-1  ^ 


(2) 


S.-  nA 


•  I  (a^  “  A) 
1-1  ^ 


We  will  then  let  A  become  a  random  variable  that  can  take  on 
uniformly  any  value  In  the  Interval  "  ^  i  ^  i.  •  Diff«rent 
choices  of  A  give  rise  to  different  rounding  errors  -  S  and 
the  variance  of  the  distribution  of  S.  -  S  can  be  used  to  measure 

O 

the  variability  of  the  rounding  error  due  to  the  random  selection  of 
the  origin  of  the  real  numbers  a^^  with  respect  to  that  of  the 
computer. 


The  cumulative  error  from  (1)  and  (2)  Is 

(3)  S.  -  S  -  I  [(a.  -  A)*  -  (a.  -  A)] 

^  1-1  ^  ^ 

Let  be  the  positive  fractional  part  of  a^  and  let  be  the 
largest  Integer  not  exceeding  a^,  l.e., 

(4)  a^  -  +  f 
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Denoting  by  r^  the  error  of  the  1^^  term,  we  have 

(5)  r.  -  [(a.-A)*  -  (a.-A)]  -< 

'l-(fi-A) 

If  -  y  <  A 
1 

XX  X 

-(f^-A) 

If  -j  +  f. 

To  prove  the  above,  we  note  that  f^-A  - 

(a^-A)  +0^ 

-  j  i  f  -A  ^  +  y  then  (a^-A)  la  rounded  to  o^.  Hence  a^-A  Is 
rounded  down  If  ~  ^  ^  otherwise  rounded  up. 


Denoting  expected  value  by  E,  we  have  by  direct  evaluation 


4 


(6)  E(r^)  -  /  r^dA  -  0 

"  2 


Assume  f^  f j  t  then 


ECriXj)  - 


-  ~ 

/  ^  2  2 

”2  ■ 
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-  /  (f^fj  -  A(£^  +  ^j>  +  ^ 


*  li  ^ 


(d-fj^-fj)  +  2A]  dA 


^  fl  ^ 

"  2^^1 


[-f^  +  A)  dA 


Performing  Indicated  Integration  yields: 

<7)  E(r^tj)  -  +  5I 
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which  is  one-half  the  2^^  order  Bernoulli  Polynomial  in  For 

£j<,£^  we  also  get  (7) .  Note  that  the  individual  errors  r^  and  r^ 
are  not  Independent  of  one  another. 

It  now  follows  that 

(9)  E(S)  -  S 

(10)  E(S-S)2-E(f  I  r  r  )-ij  f  [If  _f  |2_|f  £  1+  Ij 

i-l  J-1  ^  J  ^  i-1  j-1  ^  J  1  J  0 

-  2 

The  usual  value  of  variance,  E(S-S)  •  n/12 ,  will  result  if  we  further 
assume  f^  are  Independently  drawn  from  uniform  distributions  on 
[0  <  f^  <  1]. 

Theorem;  If  the  fractional  parts  of  all  a^  are  equal  to  each  other* 
then  each  term  of  (10)  is  maximum  for  0  £  f^  £  1'  and 

(11)  Max  ZiS-sf  -  i  X  !  4)  "  4  • 

2  1-1  j-1  ^ 

From  (10)  we  have  an  interesting  inequality*  namely  for  all  f^ 

(12)  v(f)  -  |  I  I  >  0 

^  i-l  J-l  ^  ^  1  J  0 

This  function  is  not  convex  even  for  n-2*  since  f^^^  -  (  j*  0)  and 
f  -  (-  7,  0)  yields  V(f°)  -  V(f^) 
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f  +  f  3 

V( - j - )  ■  V(0)  ■  Yj  •  There  appears  to  be  no  obvious  direct  way 

to  establish  that  V(f)  ^  0  for  all  0  ^  £  1.  Our  development  shows 
V(f)  to  be  a  variance  and  this, of  course, constitutes  an  Indirect  proof. 
We  can  replace  (12)  by  a  convex  realization:  Assume  f^  ^  for  all 

1,  then  the  problem  of  finding  Min  V(f)  can  be  rewritten: 

9  n2 

(13)  Find  Min  [V(f)]  -  ^  +  I2  "  +  (^-3)^2  + 

+...+  (n-2k+l)  f.+...-(n-l)f  ] 

K  n 


subject  to 

(14)  f^  >  fj  .  .  .  >  f„ 

(15)  0  £  f^  <  1 

Formally  (13),  (14),  (IS),  Is  a  poaltlve  definite  quadratic  program. 
Fortunately,  as  we  shall  see  this  can  be  solved  by  classical  calculus 
by  Ignoring  Inequalities  (14)  and  (15). 

t 

Theorem:  Equally  spaced  f^  *  (n  -  l)/n,  (1  ■  l,...,n)  yields 
Min  V(f)  ■  ^  Independent  of  n,  l.e.,  the  variance  of  the  sum  In 
this  case  Is  minimum  and  Is  the  same  as  the  variance  of  the  Individual 
terms  forming  the  sum. 

Proof:  Setting  partlala  "0  In  (13)  yields: 
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(n-1) 


(16)  < 


2(n-l)f^ 


-  2f, 


•  •  • 


2f^  +  2 (n-1) f 2 


-  2£ 


-  2f^  -  (n-3) 


-  2£, 


-  2£2  ...  2(n-l)£^_2  -  2f 


-(n-3) 


-  2£,  -  2£2 


2(n-l)£^  -  -(n-1) 


Adding  shows  the  equations  to  be  dependent.  Hence  we  may  drop  the  last 
equation  as  redundant.  Moreover,  we  can  always  translate  the  £^  so 
that  the  smallest  £^,  namely  "  0 


Rs-addlng  yields: 


2£,  +  2£,  +...2£^  ,  +  0  -  (n-1)  .  f  ■  0. 

X  X  u*x  n 

Adding  this  last  equation  to  each  o£  the  others  gives 

2n£^  -  (n  -  21  +  1)  +  (n  -  1)  -  2(n  -  1) 

(17)  £^  -  (n  -  l)/n 

Evidently  the  conditions  0  ^  £^  ^  1  and  ^  £^^2  (^^  luck) 

also  satls£led  so  that  (17)  yields  the  minimum,  namely 

_2  -  n  , 

(18)  Min  V(£)  -  "  I  (n-2j+l)£j^  -  ^  • 
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